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Electronic transport in nickel magnetic point contacts is investigated with a combination of density 
functional theory and the non-equilibrium Green functions method. In particular we address the 
possibility of huge ballistic magnetoresistance in impurity-free point contacts and the effects of 
oxygen impurities. On-site corrections over the local spin density approximation (LSDA) for the 
exchange and correlation potential, namely the LDA+17 method, are applied in order to account 
for low-coordination and strong correlations. We show that impurity-free point contacts present 
magnetoresistance never in excess of 50%. This value can raise up to about 450 % in the case of 
oxygen contamination. These results suggest that magnetoresistance in excess of 1,000 % can not 
have solely electronic origin. 

PACS numbers: 73.63.-b, 75.75.+a,72.25.-b 



I. INTRODUCTION 

Since the discovery of the giant magnetoresistance 
(GMR) in magnetic multilayers^ 2 - there has been an in- 
creasing interest in the electronic transport properties 
of magnetic materials. Besides establishing the ground 
for the manufacture of the present generation of mag- 
netic data storage devices, GMR has also provided a new 
paradigm whereby the electron spin as well as its charge 
can be used for electronic applications.— 

The GMR effect is the overall change in the electrical 
resistance of a magnetic device when an external mag- 
netic field is applied and it is associated to a change in 
the internal magnetic structure of the device. For in- 
stance, when the magnetic moments of the two magnetic 
layers forming a spin valve are aligned parallel to each 
other from their zero-field antiparallcl alignment, the re- 
sistance across the junction drops. This is the principle 
behind the present generation of read-heads for hard- 
drives. However, in order to reach storage densities of 
the order of Terabit/in 2 a substantial down-scaling of 
the read/write devices is needed. One possible avenue to 
this goal is offered by magnetic point contacts (MPCs), 
with typical cross sections approaching the atomic scale. 
In these the electron coherence length is greater than the 
size of the point contact and the electronic transport is 
ballistic. 

The experimental landscape for magnetic point con- 
tact transport is rather rich and controversial. Viret et 
al. have performed mechanically controlled break junc- 
tion experiments in cryogenic vacuum and shown that 
the ballistic magnetoresistance (BMR) in a one-atom- 
thick MPC can reach up to 40 % ,— in good agreement 
with theory^ Other experiments using slightly larger con- 
strictions and performed in air showed large ballistic 
magnetoresistance (LBMR) ratios reaching up to 300 % 
i ti7i?i? Finally claims of a huge ballistic magnetoresis- 
tance (HBMR) effect have been recently made with val- 
ues up to a few thousand percent. 10 : 11 

To date there is a long-standing debate around the 
origin of the LBMR and HBMR in MPCs, which is 



not completely resolved. On the one hand, it has been 
argued that magnetic field-induced mechanical effects 
can produce large magnetoresistance (MR). In fact, ei- 
ther magnetostriction ) 12 ' 13 dipole-dipole interactions^ 
and magnetically induced stress relief may cause the 
compression of the nano-contact once a magnetic field 
is applied. This enlarges the cross section of the MPC 
and consequently the resistance of the junction decreases. 
On the other hand, there is also evidence that mechan- 
ical effects alone are not able to account for whole MR. 
Specifically Garcia et. al. have shown that the behaviour 
of MPCs does not comply with mechanical changes, in 
particular with magnetostriction 10 . Hence the question 
on whether or not electronic only effects are sufficient to 
explain LBMR and HBMR remains. 

Recently Garcia et al.— have proposed that the pres- 
ence of impurities, in particular oxygen, might be related 
to HBMR (exceeding 1,000 %). This is an interesting hy- 
pothesis that concerns not only experiments conducted in 
air but also those in ultra high vacuum conditions, since 
the samples usually become contaminated after only a 
couple of hoursj 17 Importantly at these extremely small 
dimensions even a single impurity might have a large ef- 
fect on the current flowing through the device. In par- 
ticular, it is also well know that Ni is extremely reactive 
to oxygen and that bulk NiO is an insulator. Assuming 
that the insulating state of NiO persists in some form 
down to the atomic scale, then an oxygen contaminated 
Ni point contact can behave as an atomic scale tunnelling 
junction and perhaps sustain LBMR or even HBMR. 

Large MR in ballistic tunnelling junctions are not un- 
common in bulk junctions grown epitaxiall y 18 ' 19 and val- 
ues in excess of 300% have been demonstrated. Interest- 
ingly theoretical calculations have shown that the con- 
ductance is highly resonant in the two-dimensional Bril- 
louin zone orthogonal to the transport direction . 20 : 21 In 
particular, minority spins present high conductance in 
small regions away from the T point and low conductance 
everywhere else. One may then speculate that such reso- 
nant condition may be further exploited in quasi- ID sys- 
tem, where the transversal Brillouin zone collapses into 
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a single point. 

Numerical simulations have an important role to play 
in addressing these issues. In particular it is necessary 
to determine with a high degree of accuracy the under- 
lying electronic structure of the MPC and from that to 
calculate the electronic transport. This can be achieved 
by using density functional theory (DFT}2£ in its Kohn- 
Sham forrn^ 3 - and the non-equilibrium Green function 
(NEGF) formalism i 24 i 25 However, for 3d transition met- 
als in presence of oxygen it is crucial to choose correctly 
the exchange and correlation (XC) potential, since this 
can significantly change the transport properties of the 
devicej 26 i 27 Most notably, transition metal oxides are 
poorly described by the local spin density approxima- 
tion (LSDA) and corrections are needed. A particularly 
useful one is offered by the LDA+U scheme ; 28 ' 29 where 
the LSDA functional is replaced by an Hubbard-like en- 
ergy for those atomic orbitals for which strong correlation 
cannot be neglected (the 3d shell in this case). 

This paper is organised as follows. In the next sec- 
tion we describe the method used to calculate the mag- 
netoresistance as well as the two point contact arrange- 
ments investigated in our calculations. In section ITTT1 we 
present our results for the electronic transport proper- 
ties of impurity-free MPCs using the LSDA and LDA+/7 
functionals. In section HVl we discuss the role of oxygen 
impurities and finally we draw our conclusions. 



II. METHOD 

Transport calculations are performed with 
Smeago P ! 31 our electronic transport code, which 
uses a combination of DFT— iS 3 - and NEGF method 2 ^ 5 - 
to accurate predict the I-V characteristics of nano-scale 
devices. Smeagol is based on the DFT code SIESTA, 32 
which provides the Kohn-Sham Hamiltonian over a basis 
set of numerical atomic orbitals. 

We ideally divide the device under investigation into 
three distinct regions^ namely a left and a right semi- 
infinite electrode and a central scattering region where 
the potential drops^ In the case of MPCs the scattering 
region consists of the atoms pictured in either FigQ] or 
Figj^l The main quantity in our transport calculation 
is the retarded Green function of the central scattering 
region evaluated at an energy E, 

G° (E) = lim [(E + i v ) S -H°- (E) - (E)] , 

(1) 

where H a is the Kohn-Sham Hamiltonian^ 3 - for spin a 
(a =|, |), S is the overlap matrix, (^r) are the self- 
energys for the left and right lead respectively. The self 
energy contains information about the electronic struc- 
ture of the semi-infinite electrode as well as the coupling 
to the scattering region . 25 ' 30 ' 31 ' Here we consider the 
two-spin fluid approximation^ 5 - whereby the two spin 
components of the current add in parallel. In this par- 
ticular case this approach is justified by calculations per- 



formed previously by Ying et al. r& who successfully mod- 
eled a constricted junction using a mean-field Heisenberg 
model. 

A self-consistent procedure based on the non- 
equilibrium Green function has been develope d 30 ' 34 for 
calculating both G(E) and the charge density p. Once 
convergence has been achieved the energy- and bias- 
dependent transmission coefficients for each spin com- 
ponent can be calculated by using 



T(E,V) a = Tt[T l GT r G^] 



(2) 



with r L/R (E) 



J L/R 



L/R 



Finally, the spin po- 



larised conductance is simply proportional to T(E) eval- 
uated at the Fermi level Ep 



G = - [TT (E f ) + {Ef 



(3) 



We model the magnetoresistance in the MPC using 
the typical spin-valve scheme, i.e. we assume that in the 
absence of an external magnetic field, the magnetisation 
vectors of the two leads are opposite to each other in an 
antiparallel alignment (AA). Therefore, for zero field a 
sharp domain wall is formed inside the MPCj 3 -^ When a 
magnetic field is applied, the magnetic moments of the 
leads align in the field parallel to each other (PA) and 
the domain wall is eliminated from the junction. The 
low-bias "optimistic" magnetoresistance ratio -Rgmr can 
then be calculated by using the conductances of the AA 
(G AA ) and the PA (G PA ) 



R 



GMR 



G PA — G 
G AA 



AA 



(4) 



Our simulations are performed for atomic size point 
contacts, this arrangement was initially proposed by 
Viret et al..— The structure of a Ni break junction close 
to the rupture point is modelled by two Ni pyramids ori- 
ented along the [001] direction. These are formed from 
fee bulk Ni, with a lattice constant equal to 3.46A and 
they sandwich one extra nickel atom in such a way that 
a three atom-long Ni chain bridges the two Ni surfaces 
(see figure [T]) . The whole simulation cell consists of 56 
nickel atoms and the basis set comprises of double-^ (DZ) 
orbitals for the 4s, Ap and 3d orbitals, an additional po- 
larisation orbital is included for the 4s ^ The real space 
grid is given by an equivalent plane wave cut-off of 500 Ry 
and we use 64 complex energy points for integrating the 
charge density^ 30 Finally we use periodic boundary con- 
ditions in the transverse direction and we sample over 8 
/c-points in the 2D Brillouin zone. 

Before the transport calculations are carried out we 
perform conjugate gradient structure relaxation of the 
MPCs geometry. The two left-most and the two right- 
most atomic planes in the unit cell are kept fixed at the 
Ni fee bulk positions, while the middle atoms are free to 
relax. Once the relaxation of a particular arrangement is 
completed we then increase the separation d between the 
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Figure 1: (Color online) Ball and stick diagram of a nickel 
point contact formed by two pyramidal tips joined by a single 
Ni atom. The eleven atoms included in the red box corre- 
spond to the tip atoms (three atoms forming the atomic chain 
and four atoms on either side forming the apexes), d is the 
distance between the two outermost planes. 



two outermost planes (see figure [TJ and another relax- 
ation is performed. Thus we find the energy minimum as 
a function of d. The final distance between the plane of 
atoms forming the apex and the first atom in each side of 
the chain is 1.60 A and the inter-atomic distance within 
the chain itself is 2.24 A. Finally the distance between 
the plane of the apex and the first neighbouring plane of 
bulk atoms is 1.76 A. 




III. IMPURITY-FREE MPCs 
A. LSDA calculations 

We start our analysis by looking at impurity-free 
MPCs. The projected density of states (PDOS) can pro- 
vide some initial insight into the character of the elec- 
tronic states lying close to the Fermi level. In FigJ3] we 
present the density of states projected onto the s and d 
orbitals of the entire MPC and on those of the atoms in 
the tip. 
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Figure 2: (Color online) Ball and stick representation of the 
relaxed atomic arrangement of a nickel point contact with one 
oxygen atom brigding the gap between the two tips. Color 
code: grey (blue) Ni; black (red) O. 



Figure 3: (Color online) PDOS for impurity-free MPC: a) en- 
tire simulation cell, b) 11 atoms in the central region (apexes 
plus atomic chain) , c) apex atoms, d) the 3 Ni atoms forming 
the monoatomic chain. 



We also consider the possibility of oxygen impurities 
in the MPC, by replacing the middle nickel atom of the 
single-atom chain by an oxygen atom (figure [J) . In this 
case we perform structural relaxation to find the most 
energetically favourable arrangement of the device. For 
this calculation we use, in addition to the basis set of 
Ni, a double-C polarised basis for s and p orbitals of O. 
Atomic relaxation is performed from two different initial 
geometries. The first has the oxygen atom placed along 
the axis of the two pyramids, while in the second this is 
displaced perpendicularly to the axis as a continuation 
of the fee lattice. We find the same energy minimum 
for both the arrangements with the O atom in a straight 
configuration. In the transport calculations the relaxed 
coordinates are always considered. 



First we note that the DOS of the entire MPC is sim- 
ilar to the one of the apexes and resembles that of bulk 
nickel^ In contrast the PDOS for the three-atom chain 
(Fig. [3ji) shows the effects of low-coordination, since it is 
much sharper than that of the rest of the MPC. We then 
expect peaks in the transmission coefficients. These are 
presented in Figure 2] for both the magnetic configura- 
tions of the MPC. In the PA configuration we observe 
minority-dominated transmission with T^(E-p) w 2.5. 
The majority spins contribute with about (E-p) ps 1 
to give a total conductance at zero bias around 3.5 e 2 /h. 
In the AA, albeit still dominant, the minority peak is 
largely suppressed and T(E-p) takes value similar to that 
of the majority spin in the PA. In fact, we observe a 
plateau at T (E) w 1 for E > Ep and both the spins. 
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This originates from the presence of nearly un-polarised 
s electrons whose transmission is rather insensitive to the 
magnetic configuration of the MPC. The resulting MR 
ratio is approximately 20%. 
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Figure 4: (Color online) Transmission coefficients as a func- 
tion of energy for Ni point contacts formed by a single atom 
chain: a) parallel and b) anti-parallel magnetic configuration. 

The presence of unpolarised s electrons sets a serious 
limitation to the maximum MR achievable. As we have 
just shown s electrons contribute with about 2e 2 /h to 
the conductance of both the AA and the PA. The re- 
maining contribution originates from d electrons. This 
is sensitive to the magnetic configuration, but it cannot 
exceed 5e 2 /h for each spin, since at most 5 d-bands can 
cross Ef. This gives us an upper bound for the MR of 
about 250%. In practice however, d electrons are strongly 
backscattered and they never contribute with more than 
1.5 to the transmission coefficient. This indicates that 
huge MR is hardly possible in impurity-free MPC. 



B. LDA+lf Calculations 

Although useful and accurate in many situations, the 
LSDA severely underperforms for a variety of systems 
where strong electron correlations are important. Bulk 
Ni does not figure in this class of materials and in fact 
its electronic properties are extremely well described by 
the LSDAi2£ However the same cannot be said of Ni in 
atomic-size chains, where the localisation of the 3d or- 
bitals can be enhanced by the low coordination. 

In figure [5] we present a schematic diagram of the 3d 
orbitals of a transition metal mono-atomic wire aligned 
along the z direction. While the d z 2 , d yz and d xz orbitals 
are parallel to the chain and overlap considerably (note 
that in this arrangement symmetry forbids s-d hybridi- 
sation for all the 3d orbitals except for d z i). The d xy 
and d x 2_ y 2 are perpendicular to the chain and remain 
localised. The key point is that, in contrast to noble 
metals i 17 i 30 the localised d orbitals in transition metals 




Figure 5: (Color online) Schematic representation of the 
atomic orbitals responsible for the electronic structure of a Ni 
mono-atomic wire. The d x 2_ y 2 and d xy orbitals are perpen- 
dicular to the axis of the wire and therefore weakly coupled. 



appear close to Ep and may affect drastically the low- 
bias electron transport. For this reason it is crucial to 
describe their position correctly, and therefore we need 
to consider XC potentials beyond LSDA. 

Unfortunately, most of the corrections tend to deteri- 
orate the LSDA electronic structure of bulk Ni. In the 
case of point contacts, the problem then becomes that of 
finding a good functional for the low connectivity apex 
region which at the same time does not alter the elec- 
tronic structure of the planes at the edges of the cell. 
The LDA+U metho d 28 ' 29 is particularly suited for this 
purpose. The approximation consists in replacing the 
LSDA XC energy with an on-site Hubbard-type energy, 
which applies only to those atomic shells that need to be 
corrected (the d shell in this case). The LDA+C/ poten- 
tial depends on two parameters, the Coulomb repulsion 
U and the exchange J. These are associated to their cor- 
responding atomic counterparts, although in a solid the 
atomic values might be severely corrected by local screen- 
ing. We can thus consider U and J non-vanishing only in 
the region of the point contact where the Ni atoms have 
low coordination (the apexes and the tip). 

In order to understand the effect of the LDA+C/ on 
the electronic structure of the MPCs and to extract re- 
alistic values for U and J, we have first performed band- 
structure calculations for an infinite one-dimensional 
nickel chain (lattice spacing 2.24 A). For these calcu- 
lations we have used the LSDA, LDA+U with various 
choices of U and J and the atomic self-interaction cor- 
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rected (ASIC) LSDA method^ ASIC is a fully ab ini- 
tio scheme, which corrects for the atomic part of the 
LSDA self-interaction and it is proved to perform ex- 
tremely well for a broad range of materials including 
transition metal monoxides. Recently it was also used 




Figure 6: (Color online) Band structures of a one-dimensional 
nickel chain (lattice spacing 2.24 A), a) LSDA, b) ASIC and 
c) LDA+U with [7=6 eV and J=l eV. The panels to the left 
(right) correspond to majority (minority) spins. 

for transport calculation s 26 ! 27 with considerable improve- 
ments over LSDA. The bandstructure around the Fermi 
energy thus obtained are presented in figure [5] 

Generally speaking the 3d manifold splits into three 
distinct bands, a broad hybrid s-d z 2 band approximately 
8 eV wide and two doubly degenerates bands respectively 
of d yz , d xz and d xy , d x 2_ y 2 symmetry. The first of these 
degenerate bands is about 2.5 eV wide (it is located at 
energies between -3 eV and -0.5 eV for the LSDA major- 
ity spin band), which is a direct consequence of the fact 
that the d yz and d xz orbitals are aligned along the chain. 
In contrast the d xy , d x 2_ y 2 band is only 0.5 eV wide (be- 
tween -1.5 eV and -1 eV for the majority spin of the 
LSDA bands). The spin-split is of about 1 eV similarly 
to bulk Ni. The main effect introduced by the LDA+C7 
and the ASIC is that the d xy , d x 2_ y 2 band are shifted to- 
wards lower energies, thus moving them away from Ep. 
Taking the ASIC bands as reference we have fixed the 
Coulomb and exchange LDA+C/ parameters respectively 
to C/=6 eV and J=l eV. Importantly we note that the 
s-d z 2 band anti-crosses at around its band-centre and it 
is always present at the Fermi level regardless of the spin 
orientation and the XC functional used. One can thus 
expects a large portion of the conductance through the 
chain to be insensitive to the magnetic state. 

We then move to the calculation of the transport prop- 
erties of the MPC of Fig. [TJ Here we apply the LDA+C7 
corrections only to the three atoms forming the mono- 
atomic chain. Test calculations where the corrections are 
extended to the apexes yield essentially the same results. 
The zero-bias transmission coefficients for the PA and AA 



and for the different spins are presented in figure [3 little 
difference is observed when the transmission coefficients 
are calculated with LSDA (see FigfJ|. In the PA case 
the transmission for the majority spins are close to unity 
for a wide range of energies, indicating less hybridization 
between s and d orbitals, while the minority spin still has 
a peak of T w 2.5 at Ep. In general the main difference 
with the LSDA results is that now the transmission for 
energy below Ep is reduced and the peaks are broadened. 
This is the result of the downshift of the d orbitals due 
to the on-site U correction. 

Also the transmission for the AA is similar to that cal- 
culated with LSDA, except for a general reduction of T. 
This is now rather close to unity at the Fermi energy for 
both the spin configurations (note that since the domain 
wall is located between the second and third atom in the 
chain, the two spin directions are not degenerate for the 
AA). 
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Figure 7: (Color online) Transmission coefficients as a func- 
tion of energy using LDA+t/ for Ni MPC. a) Parallel and b) 
anti-parallel alignments of the magnetic moments of the elec- 
trodes. The values of U and J are 6 eV and 1 eV respectively. 

The results of figure [7] can be understood by noting 
that the MPC investigated is essentially formed by a lin- 
ear mono-atomic Ni chain sandwiched between two Ni 
surfaces. The transmission coefficient thus reflects closely 
the bandstructure of the chain itself. Consider first the 
majority spins for the PA. The large plateau in (E) for 
E > Ep is associated with the broad s-d z 2 hybrid band. 
This is the only one present in the bandstructure of the 
Ni chain in that energy range. Below Ep the contribu- 
tion of the remaining d orbitals becomes important and 
T exceeds unity. Such a contribution is present at Ep 
for the minority spins, for which 4 bands cross the Fermi 
level. One thus expects the conductance of the minority 
spins to reach up to 4e 2 /h. However scattering with the 
pyramid-shape electrodes results in a value of only about 
2.5e 2 /h in the actual MPC. 

Similarly the transmission for the antiparallel configu- 
ration can also be related to the bands of Fig. [5] In this 
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case the MPC can be thought of as two Ni chains with 
opposite magnetization direction, coupled to each other 
through an atomically sharp domain wall. Therefore 
electrons propagating in the majority (minority) band 
across the first half of the chain, must propagate as mi- 
nority (majority) in the other halfi^ Consequently the 
transmission is large at those energies where bands with 
the same orbital character appear for both spin direc- 
tions. The only band that satisfies this criterion for en- 
ergies around Ep is the hybrid s-d z 2 , which explains why 
T(Ep) ~ 1 for the AA. In contrast the remaining "pure" 
d bands are narrow and exchange split and contribute 
little to the total transmission. 

The GMR in this case is enhanced with respect to 
that calculated with LSD A and can reach up 60 %. This 
is in agreement with calculations performed with differ- 
ent ab initio methods^^ and experiments for atomically 
thin junctions. - However, this value is still far from those 
observed in some of the experiments where the HBMR 
is measured . 10 ' 11 The broad s-d z 2 hybrid band is both 
highly conductive and crosses the Fermi energy for the 
two spin sub-bands. This makes it difficult to compare 
our ab initio results with that class of experiments. More- 
over it is also important to note that the dominant s 
orbital character of the s-d z 2 band makes such conduc- 
tance channel rather robust against conformal changes 
of the point contact. In fact calculations with different 
MPC geometries give rather similar results. Thus, one 
must conclude that impurity free MPCs can not produce 
HBMR. 

Next we consider the case of MPCs contaminated with 
strongly electronegative impurities such as oxygen. Our 
expectation is then to remove the contribution to the 
conductance originating from the s orbitals. 
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Figure 8: (Color online) Band structures (left panel) and 
projected density of states (right panel) for an infinite one- 
dimensional NiO wire in the antiferromagnetic configuration 
calculated with the LSDA. In this case the two spin sub-bands 
are degenerate. 




Figure 9: (Color online) Band structures and projected den- 
IV. OXYGENATED Ni MPCs s jty of states for an infinite one-dimensional NiO wire in the 

ferromagnetic configuration calculated with the LSDA. The 
A. LSDA Calculations left ( ri ght) panels are for the minority (majority) spins. 



Following the suggestions of Garcia et alJ£ we consider 
the effects of oxygen contamination in MPCs. We start 
by calculating the LSDA bandstructure of a mono-atomic 
NiO linear chain. In doing this we follow the same idea 
of the previous section, i.e. that the transport properties 
of the MPC can be understood from the electronic struc- 
ture of the one dimensional chain associated to the three 
atoms forming the narrower part of the constriction. 

In this case the chain associated to the MPC of fig- 
ure is a NiO chain with a bond length of 1.8 A. The 
bandstructure and the corresponding PDOS for the an- 
tiferromagnetic and ferromagnetic ground state are pre- 
sented respectively in figures [5] and [5] The LSDA total 
energies of these two configurations are identical, which 
means that the chain is actually paramagnetic. Metal- 
licity is found regardless of the magnetic state and it is 
due to a hybrid s-d z 2 band. This band receives contribu- 
tions from the O p z orbitals and has a spin-split of about 



0.5 eV in the ferromagnetic configuration with the ma- 
jority component completely filled. Thus, if one assumes 
that the transport through the chain is completely dom- 
inated by such a highly conducting band, then the ferro- 
magnetic one-dimensional NiO will effectively behave as 
a half-metal with a potentially large GMR. However in 
the antiferromagnetic configuration such band splits into 
two narrow bands respectively at the Fermi level and at 
about 0.6 eV below Ep, precluding the possibility of a 
strong suppression of the current for the AA. 

With these results in hand we move on to calculating 
T(E) for the structure shown in figureEJ these results are 
shown in figure [TU1 The transmission coefficients at the 
Ep for the PA configuration are large for the minority 
and small for the majority spins (respectively 0.5 and 
2.5), this is consistent with our discussion based on the 
NiO chain. However T(Ep) is rather large for both spins 
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Figure 10: (Color online) Zero-bias transmission coefficients 
of oxygen contaminated Ni point contacts calculated with the 
LSDA: a) parallel and b) anti-parallel alignments. In the 
PA the solid (dashed) line represents the majority (minority) 
spins, while the two spins are degenerate in the A A. 
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Figure 11: (Color online) Band structures and projected den- 
sity of states for an infinite one-dimensional NiO wire in 
the ferromagnetic configuration calculated with the LDA+C7. 
The left (right) panels are for the minority (majority) spins. 
The values of U and J are set to 6 eV and 1 eV respectively. 



in the AA because of the contribution from the Ni s-d z 2- 
O p z band. This gives an overall conductance larger than 
that of the PA and a tiny negative GMR of approximately 
-6 %. It therefore appears that oxygen impurities are 
unable to explain HBMR. 

This negative result however should be taken with 
some care. As mentioned in the introduction LSDA 
poorly reproduces the electronic structure of standard 
Mott-Hubbard insulators such as NiO. One may argue 
that the same problem appears in these low dimensional 
oxides and therefore our results need to be tested against 
the LDA+J7 scheme. 



B. LDA+U calculations 

As in the previous sections we begin with the LDA+(7 
calculation of an infinite NiO chain. Here we consider the 
same values of U and J used for the case of the Ni-only 
chains. The configuration presenting the lowest energy 
turns out to be non-magnetic, 43 although both ferromag- 
netic and antiferromagnetic solutions can be stabilized. 
Interestingly several antiferromagnetic states, both con- 
ducting and insulating, have been obtained with the con- 
vergence being extremely sensitive to the initial orbital 
occupation. 

In figure QT] we present the bandstructure for the ferro- 
magnetic configuration. These appear significantly differ- 
ent from their LSDA counterpart, since NiO is now a half- 
metal. Moreover the s electrons are no longer present at 
the Fermi level, which in contrast is characterized by or- 
bitals orthogonal to the chain axis, namely the oxygen 
p x and p y and the Ni d xz and d yz . 

We then proceed with calculating the transport 
through the point contact. When the structure of figure 



[2] is considered the non-magnetic solution is no longer 
found. The transmission coefficients calculated using 
LDA+J7 are shown in figure IT2l The PA configuration is 
qualitatively similar to the LSDA result (see figure [TT)a) 
however the transmission for minority spins is slightly 
higher while that for majority is lower (T ~ 0.2). The 
majority d states are lower than in the LSDA case and al- 
though the transmission coefficients at E-p are non-zero it 
does resemble the half-metal behaviour observed in the 
infinite chain. The residual conductance at E-p is then 
attributed to direct tunnelling at the apexes, which is 
substantial for these small separations. 

For the AA configuration the LDA+C/ results are in 
stark contrast to the LSDA ones and in particular there 
is a drastic reduction of the transmission at Ep. This 
is somehow expected since the contribution of the s elec- 
trons to the transmission has been severely reduced. The 
resulting magnetoresistance reaches a much higher value 
of 450 % in good agreement with results obtained with 
by using hybrid XC functionals. 42 



V. CONCLUSION 

In summary, we have investigated the possibility of 
LBMR and HBMR in Ni point contacts using a combi- 
nation of DFT and NEGF method. In particular we have 
explored the dependence of the MR on the XC function- 
als used and the effects of oxygen contamination. For 
defect-free point contacts both LSDA and and LDA+i7 
calculations agree in setting an upper bound to the MR in 
the range of 50 %. This rather small value is essentially 
due to the high transmission associated to an unpolar- 
ized s-d z 2 band, which is always present at the Fermi 
level. This is rather insensitive to conformal changes of 
the MPC and contributes 2 e 2 /h to the conductance re- 
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Figure 12: (Color online) Zero-bias transmission coefficients 
of oxygen contaminated Ni point contacts calculated with the 
LDA+f/ : a) parallel and b) anti-parallel alignments. In the 
PA the solid (dashed) line represents the majority (minority) 
spins, while the two spins are degenerate in the AA. 

gardless of the magnetic state of the device. 

In contrast, oxygenation can drastically change the pic- 
ture. In particular O p states can hybridize with the s-d z 2 
band shifting it away from Ep. However, in this case the 
results are strongly dependent on the exchange and cor- 
relation functional, with LSDA giving a small negative 



MR while LDA+U gives a positive MR of about 450 %. 
Since the well-known problem of LSDA of dealing with 
transition metal monoxides, we believe that the LDA+Z7 
results are the most relevant. Importantly even in this 
case the MR is still much smaller than those measured 
in experiments showing HBMR. The main point is that, 
although strongly suppressed, there is still substantial 
current for the AA, attributed to direct tunnelling of the 
apexes. Certainly larger oxygenation, or the formation of 
long NiO chains may suppress this leakage current, how- 
ever we believe that highly oxygenated configurations are 
rather unlikely. This leads us to conclude that HBMR 
may have additional origins besides electronic considera- 
tions. 
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